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Abstract 

We derive relations between spin-dependent nuclear and nucleon gi and g2 structure functions 
within the nuclear impulse approximation, which are valid at all Q^, and in both the resonance and 
deep inelastic regions. We apply the formalism to the specific case of the deuteron, which is often 
used as a source of neutron structure information, and compare the size of the nuclear corrections 
calculated using exact kinematics and using approximations applicable at large Q^. 
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1. Introduction 



The study of nuclear effects in deep inelastic structure functions has by now a long 
and rich history. The importance of nuclear structure in high-energy scattering of leptons 
from nuclei was most prominently thrust into the limelight by the "nuclear EMC effect" 
discovered by the European Muon Collaboration (EMC) which found much larger than 
anticipated differences between structure functions of heavy nuclei and those of deuterium. 
Many theoretical and experimental studies of nuclear effects on structure functions followed, 
and over the years an extensive phenomenology has been developed, even though many 
questions about the origin of the effect still remain (for reviews see e.g. Ref. j2|). 

While the early studies of nuclear medium modifications focussed on heavy nuclei, where 
the magnitude of the effects is largest, it has since been realized that resolving the dynamical 
origin of the EMC effect requires understanding of light nuclei as well. Until recently, an 
anomalous situation existed whereby nuclei with A ^ 4, for which the most detailed micro- 
scopic theoretical calculations were possible, had the least empirical information available. 
A dedicated experiment at Jefferson Lab j^, with the goal of making precise measurements 
of the nuclear dependence of structure functions in a variety of light nuclei, has recently 
completed taking data which are currently being analyzed. 

Light nuclei, such as deuterium or ^He, are also often used as effective neutron targets in 
experiments seeking to extract information on the structure of free neutrons. This has been 
especially relevant for neutron spin structure studies, for which polarized '^He is commonly 
used. If one is to extract reliable information on the structure of the neutron, it is important 
that the nuclear corrections are properly accounted for, especially given the ever increasing 
precision of modern experiments. 

In recent years several theoretical analyses have been devoted to the nuclear EMC effect 
m A = 2 and A = 3 nuclei, both for unpolarized [4-16] and polarized [17-25] structure 
functions, which have quantified the effects of binding, Fermi motion, as well as relativity 
and nucleon off-shell effects. Most of these studies have focussed on nuclear effects in the 
deep inelastic scattering (DIS) region, where the exchanged photon's virtuality and the 
mass W of the final hadronic state are both large (Q^, ^ 1 GeV^). 

On the other hand, considerable data have been accumulated recently for W ^ 2 GeV, 
in the region dominated by nucleon resonances. The resonance region has received renewed 
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interest partly due to the remarkable phenomenon of Bloom-Gilman duality 26|], in which 
averages of resonance structure functions have been found to approximately equal the scal- 
ing region functions, measured at much larger W (see Ref. 27| and references therein). 
Extracting information on the neutron in the resonance region from nuclear data is partic- 
ularly challenging. 

The resonance region has considerably richer structure because of the prominence of 
specific resonances in the inclusive spectrum, such as the A(1232) or the Sn A^*(1535) reso- 
nances. As well as exhibiting pronounced peaks, these resonance structures are also strongly 

dependent. At high the role of resonances is restricted by kinematics to the x ~ 1 
domain, however, at low (~ 1 GeV^) resonances actually dominate the cross section. 
The effects of nucleon Fermi motion are expected to wash out much of these structures in a 
nucleus, leading to dramatic differences between structure functions of nucleons and nuclei, 
and hence much more interesting EMC effects than in the deep inelastic region. 

The theoretical tools for the study of nuclear effects in structure functions have, for his- 
torical reasons, been developed specifically for the DIS region, usually assuming the Bjorken 
limit (where both and the energy transfer to the target — oo), in which the kinematics 
greatly simplify. To date, however, only approximate methods have been used to describe 
nuclear corrections in the resonance and low-Q^ regions, mostly using either effective polar- 
izations, or convolution approximations which are valid strictly only in the Bjorken limit. 

In this paper we rectify this problem. Working within the nuclear impulse approximation, 
in which the virtual photon interacts with a single nucleon inside the nucleus, we derive a set 
of formulas for the spin-dependent nuclear gi and g2 structure functions, which are valid over 
a broad range of kinematics. In fact, since the derivation involves an exact treatment of finite 

kinematics, and does not depend on the twist expansion, our results are valid for any Q^, 
and in both the deep inelastic and resonance regions (and even for real photoproduction). 

The general derivation is made within a relativistic framework; however, for practical 
applications we specialize to the weak binding limit (WBL), in which the nuclear matrix 
elements are systematically expanded in powers of nucleon momentum p/M, where M is 
the nucleon mass. In this limit, we find that the nuclear gf^ and (7^ structure functions can 
be written as generalized convolution of the nucleon gl2 structure functions (r = p,n) and 
spin-dependent nucleon energy-momentum distributions in nuclei. At finite Q^, the g^2 
functions receive contributions from both gl and g2, in contrast to the Bjorken limit results 
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which are diagonal for gi [19|, l20|, [24 1 . 



Our formal results are quite general and applicable to all nuclei. However, in view of the 
current importance of understanding of nuclear corrections in light nuclei, we demonstrate 
our formalism by applying it to the simplest nucleus, namely the deuteron. In Sec. 2 we 
outline the derivation of the nuclear structure functions in the weak binding limit, and 
present for the first time the complete set of formulas for polarized nuclear structure functions 
in terms of those of bound nucleons, valid at arbitrary Q^. At finite we find a breakdown 
of the simple convolution expressions for the nuclear structure functions, in which all of the 

dependence is absorbed into the nucleon structure functions, with the nucleon momentum 
distributions being functions of the nucleon momentum and energy only. In the generalized 
convolution that we derive at finite Q^, the nucleon momentum distributions depend in 
addition on the parameter 7 = \q\/qo, where q and go are momentum and energy transfer, 
respectively. 

In Sec. 3 we apply the formalism to the specific case of the deuteron, and study the 
dependence of the finite-Q^ nucleon momentum distributions on the parameter 7. The 
numerical results for the deuteron structure functions are presented in Sec. 4, where we 
focus in particular on the EMC effect in the nucleon resonance region, and contrast this 
with the effect for leading twist structure functions. Finally in Sec. 5 we summarize our 
results and discuss future work. 



2. Nuclear structure functions in the weak binding limit 

We begin our discussion of the nuclear structure function by summarizing the results 
within the framework of the relativistic nuclear impulse approximation, in which the scat- 
tering from the nucleus proceeds via the scattering from its individual nucleon constituents. 
Corrections to the impulse approximation, in the form of multiple scattering from nucleons, 
or meson-exchange currents, are typically confined to the small-a; region, and can be safely 
neglected by restricting the analysis to x ^ 0.2. Possible corrections due to the final state 
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interaction of the spectator nucleon with the hadronic debris are also not considered here. 



2 a. Hadronic tensor 

In the impulse approximation the hadronic tensor for a nucleus with four-momentum 
and spin S can be expressed in terms of the nucleon propagator in a nucleus and the virtual 
photon Compton amplitude for a bound nucleon (for more details see e.g. Refs. |8|,|9|,Q, 
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W^,{Pa, q,S) = I [dp] Tr [A^{p, Pa, S) W;,{p, q)] , (1) 

where the sum is taken over bound protons and neutrons (r = p,n), and the integration 
is performed over the nucleon four-momentum p, for which we use the contracted notation 
[dp] = d^p/{2'ji)^. In Eq. ([T]), A'^{p, Pa, S) is the imaginary part of the proton (r = p) or 
neutron (r = n) propagator in a nucleus A with momentum Pa and spin S: 

Al^ip, Pa, S) = j dt dV e'(^'ot-p-) (p^, ^| Wp{t, r) vl/;(0) \Pa, S) , (2) 

with r) the (relativistic) nucleon field operator, and a and [5 Dirac spinor indeces. The 
bound (off-shell) nucleon electromagnetic tensor VV^^(p, g) in Eq. ([1]) is a matrix in Dirac 
space, and the trace "Tr" is taken in the nucleon Dirac space. 

The analysis of Eq. ([T]) in the fully relativistic case requires the solution to the nuclear 
bound state problem [in particular the calculation of A^ip, Pa, S)], a task yet to be fully 
completed.^ The presence of nucleon spin introduces some complications in Eq. ([T]), such as 
additional Lorentz-Dirac structures (structure functions) in the hadronic tensor VV^j,(p, g) 
for the bound nucleon j^, 9|. Nevertheless, the analysis can be significantly simplified in 
the target rest frame within a nonrelativistic approximation, or weak binding limit (WBL), 
assuming typical nucleon momenta and energies to be small compared to the nucleon mass. 

Starting from the general expression for Wj^j^{p, q) for the off-shell nucleon, and performing 
a systematic expansion of the matrix elements in Eq. ([1]) in terms of \p\/M and e/M up to 
order p^/ M^ and e/M, where e = p^ — M the nucleon separation energy, one can show 
' igI . 22] that the nuclear hadronic tensor can be written in terms of the (nonrelativistic) 







^ We refer in this context to calculations based on Bethc-Salpeter and light-cone approaches [l2l. |28| 
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nucleon spectral function in the nucleus, p, S), and the bound nucleon hadronic tensor, 



Ma 



[dp] 



ti [V^{e,p,S)WUp,q)] 



(3) 



M + e 

Here Ma is the nuclear mass, and the trace "tr" is taken over nucleon polarization states. 

The spin-dependent part of the nuclear hadronic tensor is related to the antisymmetric 
part of the off-shell nucleon tensor W^^,, which depends on the nucleon polarization. In the 
WBL this component of the hadronic tensor of the nonrelativistic bound nucleon is given 
by: 



M 
p-q 



S-q 
p-q 



(4) 



where g'[2 are the spin structure functions of the off-shell proton or neutron (r = n) 
with four-momentum p = (M + e^p). The operator S in Eq. (jl]) has a structure similar to 
that of the spin four-vector (0, cr) boosted to a reference frame in which the nucleon has 
(nonrelativistic) momentum p: 



S 



(7 ■ p 

M ' 



P\^P) 

2M2 



(5) 



In Eq. ([3]) the nuclear hadronic tensor factorizes into high-energy (VV^,y) and low- 
energy (P) domains in the nuclear scattering amplitudes. The presence of spin, however, 
implies that this factorization does not in general carry over to a corresponding factorization 

22| , unless further assumptions or simplifications are 
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of nuclear structure functions 
made. 

The relation between the nuclear and the nuc 
■^He and for the deuteron was discussed in Refs. 19 



eon spin structure functions for polarized 
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23 



24( 1 . Note that the implementation 



of the impulse approximation is different in Ref. [19| and in the present approach (see also 
Ref. 22[). Here we begin from a relativistic impulse approximation, Eq. ([T]), which fully 
takes into account the nucleon spin degrees of freedom, and work within the WBL to obtain 
the nuclear hadronic tensor in terms of the nuclear spectral function, Eq. ([3]). The latter 
equation is similar to the corresponding relation in Ref. lol, but is not identical due to a 



different treatment of the reaction kinematics. In particular, in Ref. 



191 ] the struck nucleon 



is assumed to be on the mass shell, while here we take the nucleon to be off-mass-shell. 
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Further details of the derivation of Eq. ([3]) will be given elsewhere 
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2 b. Spectral function 

The nuclear spectral function V is defined similarly to the nucleon propagator in the 
nucleus A in Eq. ([2]), but involves the correlator of nonrelativistic two-component nucleon 
operators: 

P;,,(e,p,5) = I dtdVdVe^(^*-P(^'-^)) {rJir',t) C(^,0)> , (6) 

where the nucleon operator ip'^ describes the nucleon with polarization a and isospin r, and 
the average is taken over the nuclear ground state. The nuclear spin is described in Eq. ([6]) 
by the axial three- vector S. For spin-1/2 targets, S is simply the target polarization vector. 
For the spin-1 case S is defined in terms of polarization vectors as S = 'ie'^*xe™, where 
m = 0, ±1 is the spin projection along the axis of quantization. 

The nuclear spectral function can in general be calculated by inserting a complete set 
of intermediate states and computing the resulting transition matrix elements between the 
ground and intermediate states. For the deuterium nucleus, the intermediate states are 
exhausted by a single proton or neutron, and the spectral function is expressed entirely in 
terms of the deuteron wave function. However, already for A = 3 nuclei the calculation of 



the spectral function is considerably more complicated [l9l. l20[|. 

The general spin structure of the spectral function can be obtained by expanding the 
matrix V^^, in terms of the Pauli spin matrices and applying constraints from parity and 
time reversal invariance. One can then write the spectral function in the general form (for 

n 

both the proton and neutron contributions) |19l |: 

V^ie,p,S) = ^if5 I + fla-S + f^T,,S,a,) , (7) 

where / is the unity matrix and Tij = % pj — | 6ij is a traceless symmetric tensor with 
Pi = Pi/\p\ for component i of the momentum (the sum over repeated indices is implied). 

Since the spectral function is hermitian, the coefficients f[ in Eq. ([71) are real functions of 
the energy e and momentum p. The function /q gives the spin averaged spectral function, 
and is normalized to the number of protons (r = p) or neutrons (r = n) in the nucleus: 

I [dp] tr [V^^-\e,p, S)] = I [dp] = Z{A-Z). (8) 



The functions f{ and /J describe the distribution of nuclear spin amongst the nucleons. The 
integrated functions and /J determine the average nucleon polarization in the target and 
the average tensor polarization, respectively: 

i^^r = J [dp] fl , (9) 
{T.. <y,Y = [dp] n , (10) 
where we take the nuclear spin vector to lie along the z-axis. 



2 c. Master formula 



Projecting from the hadronic tensor ([3]) the appropriate structure functions, we obtain 
our "master formula" for the and structure functions of the nucleus: 



(11) 



where a, 6 = 1, 2, and again summation over repeated indices is implied. In Eq. ffTTl) x' = 
Q^/2p ■ q = x/{l + {e + 'ypz)/M) is the Bjorken variable for the bound nucleon, with 



X = (5^/2Mgo, and 7 = \q\/qo = a/1 + AM'^x'^/Q'^ takes into account finite-Q^ kinematics. 
The spin-dependent nucleon momentum distribution functions Dab can be written in terms 
of the coefficient functions fi (dropping isospin labels) as: 



1^11 =/i + 



3-7^ 
672 

(3 - 7^)p^ - 1 



l)/2 



VPz 

7 



fl + 3/2 



r 



D 



12 



- 1) 



1272 
272 



(3/1 - /2) 
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/2 + — /: 

7 
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2 , l+p|(47^-3) , 5 + 18p^7^-5p2(3^27^) ^ 
I — Ji H TTT^ /2 



/^22 



472 

3p^ - 1 . ^'^^ / f , 2 
^^^2 - — + 3/2 

=/i + (3P^ - 1) /2 



,23p^-l 



1272 
2 (3/1-/2) 



127 



672 



+ 



VPz 

7 



^ + ^+^^(^-7^/2 



(l-7')/i + 

^(3-672 + 47^) -1-272 



472 



/i 



(12a) 



(12b) 
(12c) 
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5-27^(l + 3p^)+4ph' 
1272 



i^Pl - l)/2 



(12d) 



where the velocity parameter v = \p\/M. 

The derivation of Eqs. (fT2l) makes use of the fact that characteristic values of v are small, 
and one can expand the kinematical factors in a series in v keeping terms up to Olv"^). The 
functions Dab have also been averaged over the polar angle of the nucleon momentum in 
the transverse plane {px,Py) using the independence of /i,2 of the directions of the nucleon 
momentum and the independence of the Bjorken variable x' of Px and Py. 

In general the structure functions gl 2{x' , Q"^ , p"^) of the bound nucleon are functions of 



the invariant mass squared of the nucleon, since p^ ^ M 



unknown, but has been estimated in various models j^, [o], 15, 16, 22|. In the WBL, however, 
and especially for light nuclei such as deuterium, the degree to which the nucleons are off- 
mass-shell is not large, and we can assume that the off-shell functions can be approximated 
by their on-shell values, (?! 2(2;', <5^,p^) ~ fl'i 2(2^') Q^)- 

Note that while the above derivation is valid in the WBL |£:| ^ M), Eq. (fTTl) holds 
for arbitrary momentum transfer q. An interesting feature of the distribution functions Dah 
is that they depend on the momentum transfer through the dimensionless parameter 7. 
Furthermore, nuclear effects cause the off-diagonal distributions D12 and D21 to be nonzero, 
which results in mixing of the gi and (72 structure functions in the convolution integral ( ITTil . 
In leading order in v the functions D12 and D21 are driven by tensor distribution (3]3^ — 1)/2. 
In the limit of high Q^, the parameter 7 — > 1 and the distributions (fT2|) simplify considerably. 
In particular, D12 and the convolution formula for gi becomes diagonal {i.e. there are 
no contributions from g2 to g^). However, mixing in g^ persists even in this limit 22|. 
A similar mixing of the nucleon structure functions gl and the gi nuclear structure 



lis dependence is a priori 



function was observed in Refs. 



19 



20 



24| . However, the distribution functions in Eq. ([12 
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3 i 



in the higher order terms in 



are different from the corresponding results in Refs. 
V and 7^ — 1. As discussed at the end of Sec. 2a, these differences arise from the different 
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20 



2^. 



treatments of the impulse approximation here and in Refs. 

The fact that Eq. ffTTl) can be applied to both the inelastic and quasi-elastic scattering, 
at any Q^, allows us calculate nuclear structure functions in the DIS region (at large 
and W), as well as at lower Q^, where the \ow-W resonance region plays a more prominent 
role. In this context we note that the dependence of the effective nuclear distributions 
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Eq. (fT2|) on enters through the dimensionless parameter 7, a feature which allows a 
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simple parameterization of the effective distributions over the full range of kinematics 
In the next section we apply the formal results presented here to the specific case of the 
deuteron. 



3. Deuteron spin structure functions 

For the case of lepton scattering from a deuteron, the functions /j in Eq. ([71) can be 
written in terms of the deuteron wave functions as (see also Appendix B of Ref. |22|): 

/o = 47r3 {^l + ^l) 8{e -eD+ /2M) , (13a) 
A = 47r3 (V;2 - ^p|/2) 6{e - + pV2M) , (13b) 
h = 4^ ^{tPl - v^^o^2) S{e -en+ /2M) , (13c) 

where eo = —2.2 MeV is the deuteron binding energy, and ipQ and il)2 are the S- and D-state 
momentum space wave functions, respectively, normalized such that: 

/ dpp2(^o'(p)+^2(p)) = l- (14) 

Since the deuteron is an isoscalar nucleus, Eqs. f|T3l) hold for both the proton and neutron 
distributions. The average nucleon polarization and tensor polarization in the polarized 
deuteron can then be expressed as: 

(^.) = ^-\pd. (15) 
{T,ia,) = ^{Pn-V2PsD) , (16) 

where Pd = J dpp^ip2{p) is the D-state probability, and Psd = J f^V p'^''Po{p)''p2{p) is the 
S-D interference in the deuteron. For the deuteron wave function calculated from the Paris 



potential [30|] one has Pd = 5.8% and Psd = 9.4%, while for the Bonn potential [31| 
Pd = 4.3% and Psd = 10.1%. Note that in the 7 — 1 limit, the distributions Dab in 
Eqs. (1121) for the deuteron are equivalent to those in Ref. 2^, where the functions fi were 
effectively defined including a factor (1 — p^/2M^), which in our notation is now contained 
inside the distributions Dab- 

As discussed in Sec. 2 above, at finite the deuteron structure functions receive 
contributions from both the isoscalar and g2 structure functions of the nucleon (A^ = 

10 



p + n) individually. In particular, while the contribution from (7^ to gf vanishes in the 
Bjorken limit, it is non-zero at finite Q^. To illustrate the relative importance of the 
contributions to the deuteron structure functions, we can compare the individual momentum 
distributions Dab for the diagonal and off-diagonal terms, and also as a function of 7. 

In the simplest convolution model the nuclear structure functions in Eq. (ITT]) are written 
as convolutions of the nucleon structure functions and effective nucleon light- cone momentum 
distributions: 

gtix^Q') = Dabivn) 9^ {-,Q') , (17) 
Jx y \y J 

where y = {pq + ^p^jM = [1 + [e + '~fPz)/M) = x/x'. In the Bjorken limit (7 1) the 
variable y is the light-cone fraction of the deuteron carried by the interacting nucleon. The 
effective light-cone momentum distributions Dab are obtained by integrating the functions 
Dab in Eq. ([Il): 

Dab{y,l) = J [dp] Dabie,p,^) 6(^y-l- • (18) 

In Fig. [1] we show the nucleon light-cone momentum distributions Dabiy,'^) for several 
values of 7. The results for 7 = 1 correspond to the Bjorken limit distributions. Note that 
in this limit the function D12 vanishes. The diagonal functions Du and D22 are significantly 
larger than the off-diagonal functions, but decrease in magnitude for larger 7. On the other 
hand, the distribution D12 becomes larger with increasing 7, with its magnitude reaching 
~ 5% that of Dn at 7 = 2. 

In Ref. 231] the resonance region was studied using a finite-Q^ convolution formula for 
gf similar to that in Eq. 0171) . However, the non-diagonal term arising from was not 
present. We find that this term, although small, does make a non-zero contribution at finite 
values of Q^. 

In the next section we use Eq. ( ITT]) to evaluate the effects of the smearing of the nucleon 
structure functions by the nucleon momentum distributions, and compare the calculated 
deuteron structure functions with the input nucleon structure functions in the resonance 
and deep inelastic regions. 

4. Nuclear effects 

In this section we present results for the gf and structure functions of the deuteron. 
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and compare these with the free nucleon functions. We focus in particular on the resonance 
region, with W ^ 2 GeV, which has to date received httle attention. The need to understand 
nuclear effects in the deuteron gf and g2 structure functions has arisen partly in response 



to the recent hig. 



Jefferson Lab 



i-precision data on the spin dependent deuteron structure functions from 
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33 



34|. 



For the nucleon gl and structure functions, in this analysis we consider the parame- 



terizations from the MAID unitary isobar model 
al. 36j, and from EGl in CLAS at Jefferson Lab 



or (e, e'p) reactions [35|, from Simula et 



37| . These parameterizations encompass 



the resonance as well as the deep inelastic regions. In the case of the MAID model, how- 
ever, only a few selected final states are considered, so this parameterization is expected to 
underestimate the high-H^ (or low-x) region. 

In Figs. [2] and [3] we show the input xgi and xg2 structure functions, respectively, for 
an isoscalar nucleon (A^ = p + n), at a sample value, = 2 GeV^. At this the 
most prominent structure at large x in both the g^ and structure functions is the peak 
associated with the P33 A resonance, which is negative for g^ and positive for (7^. The 
three models give qualitatively similar results here, although quantitatively there are some 
differences. At lower x the second and third resonance regions are also prominent, and here 
the MAID fit is smaller in magnitude, as expected, given that it is constructed primarily to 
describe \ow-W data. 

For comparison in Fig. [2] we also show the leading twist ( "BB" ) parameterization of 
xgi from Ref. 38[, which is smooth and does not contain any resonance structure. It is 
interesting to observe that the leading twist fit appears to go through the average of the 
resonances (with the exception of the A resonance) for the MAID fit [35||, reminiscent of 



the Bloom-Gilman duality between resonance and DIS structure functions 
resonance parameterizations 



36 



27| . The other 



37j are on average larger in magnitude than the leading 



twist curve, which may reflect the presence of the nonresonant background that is included 
in these fits. 

Using the BB parameterization 38| we also calculate the leading twist Wandzura-Wilczek 



(WW) approximation to g2, g2^'^{x,Q'^) = —gi{x,Q'^) + dy gi{y,Q'^)/y, which is shown 
in Fig.O This again displays very different behavior, both in magnitude and sign, compared 
with the g2 resonance parameterizations at large x. As with the g^ comparison, the leading 
twist curve appears to average the second and third resonance regions at intermediate x as 



12 



described by the MAID 35|] fit, but is smaller in magnitude here compared with the other 



fits 



36 



37|. 



With these nucleon structure function parameterizations, we can now investigate nuclear 
effects in both the resonance and DIS regions. For the resonance region [W ^ 2 GeV), we 
use the MAID parameterization 35|] for and (7^. The deuteron xgf structure function, 
calculated using Eqs. (fTTj) . fll2al) and (112bp with the Paris deuteron wave function, is shown 
in Fig. m at = 2 GeV^ (solid curve). Here the full calculation of xgf is compared with 
that using the Bjorken limit in the momentum distributions Du and D12 (dashed), and with 
the free nucleon structure functions (dotted). Where the resonance structures are clearly 
evident in the nucleon g^ functions, they are significantly diluted in the deuteron structure 
function. At the peak of the A resonance, for example, the magnitude of the deuteron gf 
is about half that of the nucleon g^ . A similar effect is seen for the g2 structure function 
shown in Fig. [5l 

The differences between the full results for gf2 and the Bjorken limit (Q^ 00) approx- 
imation for Dab are small where the nucleon g^2 are smooth, but become more significant 
in the vicinity of the resonance peaks. At intermediate x, around the second and third reso- 
nance regions, the full results are some 10-15% smaller in magnitude than the Bjorken limit 
structure functions. At larger x the enhanced smearing is even more dramatic, with the full 
results being up to 25-30% smaller at the A resonance peak than those with the Bjorken 
limit smearing. This behavior can be understood from the nucleon light-cone distribution 
functions in Fig. [1], which decrease in magnitude with increasing 7. At fixed = 2 GeV^, 
the third resonance region at x ~ 0.5 corresponds to 7 ~ 1.5, while for the A peak at 
X ~ 0.75 one has 7^2. Therefore, at fixed Q^, larger x implies larger 7, and hence stronger 
smearing effects in both gf and g2- 

The enhanced smearing at finite-Q^ kinematics is also evident for leading twist structure 
functions, although the effects here are less dramatic. In Figs. [6] and [7] we show the deuteron 
xgf and xg2 structure functions, respectively, evaluated using the leading twist "BB" param- 
eterization of the nucleon structure functions dX = 2 GeV^. The full calculation for gf 
and the Bjorken limit approximation are similar, and both smaller in magnitude than g^ for 
X ^ 0.8. Here to a good approximation the deuteron structure functions are related to the 
nucleon structure functions by an x-independent multiplicative factor, gf2 ~ (1 — \PD)gi2- 
At larger x, one sees the onset of a classic Fermi motion effect, whereby the ratios gf 2/912 
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increase in magnitude as x — 1. In the region 0.8 ^ x ^ 1, the finite-Q^ results are again 
smaller in magnitude than those for the Bjorken limit kinematics, implying a stronger EMC 
effect at large x, for both the gf and (7^ structure functions. 

Although the shapes of the nucleon and deuteron structure functions are very different, 
especially in the resonance region, when integrated over x the differences turn out to be 
remarkably small, particularly for the lowest moment. At large the lowest moments of 
the nucleon and deuteron gi 2 structure functions, ri 2(Q^) = J dx gi^2{x,Q'^), are to a very 
good approximation related by: 

rUQ') ^ (1 - Ipd^ r^AQ") ^ m 

where for the Paris deuteron wave function the depolarization factor relating the moments 
is (1 — |-Pd) ~ 0.91. Numerically, for the leading twist BB parameterization of the nucleon 
structure functions |38|], using the Bjorken limit nucleon momentum distributions we find 
(^J)/rf = 0.91 and ^^^Vr^ = 0.91 at = 2 GeV^. Using the full expressions in 
Eqs. ([n]) and ([12]), on the other hand, we find the ratios Tf/Ff = 0.90 and T^/T^ = 0.90. 
Thus the finite-Q^ kinematics slightly reduces the magnitude of the deuteron moments, 
which is consistent with the observation of the additional suppression at finite at large 
X in Figs. E] and El 

The ratios of the deuteron to nucleon moments in the resonance region remain similar 
to the leading twist ratios, even though the shapes of the functions here are strongly 
dependent and hence infused with large higher twist contributions. Specifically, for the 
MAID parameterization of (^f^g ^^e ratios at finite-Q^ kinematics turn out to be Tf/T^ = 0.91 
and T2/T2 = 0.92, at the same = 2 GeV^. Note that the structure functions here have 
been integrated from Xmin ~ 0.4 up to Xmax = a^thr, where Xthr = Q^/(W^thr ~ + Q^) 
corresponds to the kinematical pion production threshold, Wthr = M + m,r- The results 
using Bjorken limit kinematics differ from these only in the third decimal point. 

Overall, our results suggest that for the lowest moments of the gi and g2 structure func- 
tions, the nuclear effects in the deuteron can to very good accuracy be accounted for by 
applying the depolarization correction, (1 — |-Pd)- This will not be true, however, for higher 
moments, and certainly this approximation will break down dramatically at large x, for 
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X > 0.7-0.8. 

5. Conclusions 

For most of the history of lepton-nucleus deep inelastic scattering, the discussion of 
nuclear effects on structure functions has been confined to analysis of high-energy data within 
theoretical frameworks constructed to be valid in the limit of large and W"^ M^). 
Recent high-quality data at lower and W^, especially in the transition region where 
nucleon resonances merge into the DIS continuum, have revealed a richness of phenomena 
which had not previously been appreciated because of the lack of precision in earlier data. 
An accurate description and understanding of the new data clearly demands comparable 
advances in the theoretical tools. 

Our aim in this work has been to provide the framework for analyzing modern high- 
precision data over the full range of kinematics where they are available. To this end we 
have derived relations between spin-dependent nuclear and nucleon gi and g2 structure 
functions in the weak binding limit, which are valid at all Q^, and in both the traditional 
DIS region and the poorly explored nucleon resonance region. As a consistency check, we 
verify that our results approach the previously derived convolution formulas for the nuclear 



structure functions at large 22f|. 



We apply the formalism to the specific case of the deuteron, which is often used as a 
source of neutron structure information, and compare the size of the nuclear corrections 
calculated using exact kinematics and using approximations applicable at large Q^. We find 
that significant smearing of the nucleon structure functions occurs in regions where nucleon 
resonances are prominent, with the exact results some 10-15% smaller in magnitude than 
the Bjorken limit structure functions around the second and third resonance regions. The 
smearing is enhanced at larger x, where the full results are up to 25-30% smaller at the A 
resonance peak than those with the Bjorken limit smearing. 

The enhanced smearing at finite-Q^ kinematics is also evident for leading twist structure 
functions, although the effects here are less dramatic. At intermediate x the deuteron 
structure functions are approximately given by 2 ~ (1 ~ ^^0)01 2- However, at larger x, 
0.8 ^ X ^ 1, the finite-Q^ results are again smaller in magnitude than those for the Bjorken 
limit kinematics, implying a stronger EMC effect, for both the gf and g2 structure functions. 
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The ratios of the integrals of the deuteron and nucleon structure functions are to a good 
approximation simply related by the depolarization factor, ^i^2iQ'^) /^i,2iQ'^) ~ (1 — 2^0), 
even at low {Q^ ~ 2 GeV^). This holds for both the leading twist structure functions, 
and for structure functions in the resonance region, where the shapes of the functions are 
strongly dependent, with deviations only of (9(1%). 

In closing we should mention that our calculation is by no means complete. We have 
not considered, for example, possible modification of nucleon properties (such as masses or 
widths) in the nuclear medium, which in our framework amounts to neglecting nucleon dy- 
namical off-mass-shell effects in the nucleon structure functions. Extensions in this direction 
can be carried out using existing models for the off-shell extrapolation [sl, [21I, 22, 39]; we 



have chosen not to do so in this work in order to more clearly isolate the effects associated 
with the finite-Q^ kinematics. Similarly, we have not considered here effects beyond the 
nuclear impulse approximation, such as rescattering, meson exchange currents, or final state 
interactions. Finally, our formalism can be easily applied to other nuclei, such as polarized 
^He. This will be the subject of a future publication [291]. 
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FIG. 1: Effective nucleon liglit-cone momentum distribution functions Dab{y,j) for 7=1 (Bjorken limit), 1.5 and 2. 




FIG. 2: Spin-dependent xg"^ structure function of an isoscalar nucleon, from the fits of Refs 



(solid), 
2 GeV2. 



36( 1 (dashed), 



37l | (dotted), and the leading twist fit of Ref. 



35| 



33] (dot-dashed), at 




FIG. 3: Spin-dependent xg2 structure function of an isoscalar nucleon, from the fits of Refs. 



(solid). 



36l ] (dashed). 



35| 



37l | (dotted) , and the Wandzura-Wilczek approximation to §2 using the 



leading twist gi fit of Ref. 



3i (dot-dashed), at = 2 GeV^. 
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FIG. 4: Spin-dependent xgi structure function of the deuteron evaluated at finite-Q^ (solid) and 



Bjorken limit (dashed) kinematics, compared with the nucleon (dotted) input, from Ref. 
Q2 = 2 GeV2. 
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FIG. 5: Spin-dependent xg2 structure function of the deuteron evaluated at finite-Q^ (solid) and 



Bjorken limit (dashed) kinematics, compared with the nucleon (dotted) input, from Ref. 
Q2 = 2 GeV2. 
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FIG. 7: Leading twist xg2 structure function of the deuteron evaluated at finite-Q^ (solid) and 



Bjorken limit (dashed) kinematics, compared with the nucleon (dotted) input from Ref. 
g2 = 2 GeV2. 
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